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Abstract 

The one-loop partition function for a charged self-interacting Bose gas at 
finite temperature in D-dimensional spacetime is evaluated within a path 
integral approach making use of zeta-function regularization. For D even, 
a new additional vacuum term — overlooked in all previous treatments and 
coming from the multiplicative anomaly related to functional determinants — 
is found and its dependence on the mass and chemical potential is obtained. 
The presence of the new term is shown to be crucial for having the factoriza- 
tion invariance of the regularized partition function. In the non interacting 
case, the relativistic Bose-Einstein condensation is revisited. By means of a 
suitable charge renormalization, for D = 4 the symmetry breaking phase is 
shown to be unaffected by the new term, which, however, gives actually rise 
to a non vanishing new contribution in the unbroken phase. 
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I. INTRODUCTION 



The one-loop finite temperature effective potential with a non vanishing chemical po- 
tential has been considered quite often during the last years. The case of a free relativistic 
charged bosonic field was investigated in Refs. while the self-interacting charged scalar 

field has been studied in Recently, the issue has been reconsidered and critically ana- 
lyzed in Ref . , where use has been made of zeta-function regularization in the calculations 
and a different approach for including the chemical potential has been proposed. 

We would like to recall the importance of zeta-function regularization, introduced in 
PHT(][, as a powerful tool to deal with the ambiguities (ultraviolet divergences) present in 



relativistic quantum field theories (see for example |TT|,[T2|). It permits to give a meaning 
— in the sense of analytic continuation — to the determinant of a differential operator which, 
as the product of its eigenvalues, would be formally divergent. For the sake of simplicity we 
shall here restrict ourselves to scalar fields. In the case of a neutral scalar field, we recall that 
the one-loop Euclidean partition function, regularised by means of zeta-function techniques, 
reads [0 

InZ = In det ^ = ^('(01^^,) + \mLn)\riM' , 

where C,{s\L£,) = Tr L]^^ is the zeta function related to Ld, a second order elliptic differential 
operator, C'(0|Lb) its derivative with respect to s, and is a renormalization scale mass. 
The fact is used here that the analytically continued zeta-function is generically regular at 
s = 0, and thus its derivative is well defined. 

In this paper, we shall take full advantage of a rigorous use of this regularization in 
investigating the finite temperature effects in presence of a non vanishing chemical potential. 
In order to illustrate the new issue we shall be dealing with, we start with the partition 
function for a free charged field in IR^, described by two real components 0j. The Euclidean 
action is 

^ = 1 rfx^[0,(-Az, + m2)0,] , (1) 

where 0^ = (j)k4>k is 0(2) invariant. 

As is well known, in order to study finite temperature effects, one observes that the grand 
canonical partition function may be written under the form (we shall call this method of 
including the chemical potential, as discussed in Ref. 0, method I) 

Z^(^) = Tre-'3(^-'^«), (2) 

in which H is the Hamiltonian, Q the conserved charge operator and (3 the inverse of the 
equilibrium temperature. Making use of a path integral representation and integrating over 
the momenta, one arrives at the following recipe [||,|[^^|13|: (i) compactify the imaginary 
time T in the interval [0,j3], (ii) assume a periodic boundary condition in P, and (iii) include 
the chemical potential /i by adding to the action the term 

ineeik(l)idr(j)k - ^e^V^ , (3) 
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where e is the elementary charge. As a result, the partition function reads 



l<t>{T) = <t>(T+f3) 

where D = N + 1 and A is given by 
with 

Lij = {hr + Ln — e^fi^'j 6ij, Lm = — Aat + 



(4) 



(5) 



(6) 



in which Aat is the Laplace operator on (continuous spectrum k"^) and Lr = —d"^ 
(discrete spectrum uj"^ = ); the Laplace operator on S^. Thus, one is actually dealing 
with a matrix-valued elliptic differential non self-adjoint operator acting on scalar fields in 
X R^. In this case, the partition function may be written as ||^ 



InZp^fi) = — -Indet 



A. 



ik 



M2 



-- In det 
2 



M2 M2 



where 



L± = Lt + Li^ + e jj, ± 2ejj,{L 



N) 



IS 



Another possible factorization (see, for example, |]T2[ 

K± = Ln + Lr - ± 2iefi{Lr] 



(7) 



(8) 



(9) 



Of course we have L^L^ = K^K^ and, in both cases, one is dealing with a couple of 
pseudo-differential operators (\l'DOs), and L_ being also formally self-adjoint. 

As is clear, in any case the product of two elliptic \E'DOs, say A+ and A_, appears. 
It is known that, in general, the zeta-function regularized determinants do not satisfy the 
relation det(y4i?) = det A det B. In fact, in general, there appears the so-called multiplicative 
anomaly lHJll. In terms of F{A,B) = det (AB) /{det A det B) [||], it is defined as: 



aD{A,B) = lnF{A,B) = lndet(A5) - Indet(A) - Indet(S) 



(10) 



in which the determinants of the two elliptic operators, A and B, are assumed to be defined 
(i.e. regularized) by means of the zeta function P|. Thus, the partition function, chosen a 
factorization is given by 



In Zb( A: ,A_) = — Indet = — Indet — (11) 

= + \c'm-) + [C(0|A+) + cm-)] - \an{A^, A_) . 

Here one can see the crucial role of the multiplicative anomaly: as, in the factorization, 
different operators may enter, the multiplicative anomaly is necessary in order to have the 
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same regularized partition function in both cases, namely In Z^(y4+, = In Zj3{Bj^, B^), 
which is an obvious physical requirement. (What it is very easy to see in general is that det 
det A_ ^ det B+ det B_.) 

With regards to the multiplicative anomaly, an important result is that it can be ex- 
pressed by means of the non-commutative residue associated with a classical \E'DO, known 
as the Wodzicki residue [|T6 |. 

In the limit of zero temperature and vanishing chemical potential, the two elliptic \E'DOs 
L+ and L_ reduce to Laplace-type second order differential operators with constant potential 
terms defined in H^. In this case, the multiplicative anomaly has been computed in Ref. 

One could argue that the presence of a multiplicative anomaly is strictly linked with the 
zeta-function regularization employed and, thus, that it might be an artifact of it. In the 
Appendix we will prove that the multiplicative anomaly is present indeed in a large class of 
regularizations of functional determinants appearing in the one-loop effective action. 

The contents of the paper are the following. In Sec. 2 we reconsider the free case mak- 
ing use of zeta-function regularization and we carry out a comparison of the two possible 
factorizations. In Sec. 3, the self-interacting case is presented in the one-loop approxi- 
mation. In Sec. 4 we briefly introduce the Wodzicki residue and present a proof of the 
Wodzicki formula expressing the multiplicative anomaly in terms of the corresponding non- 
commutative residue of a suitable classical ^DO. In Sec. 5, the results of Wodzicki are used 
in the computation of the multiplicative anomaly for the interacting 0(2) model. In Sec. 6, 
the Bose-Einstein condensation phenomenon is discussed in the non-interacting case. Some 
final remarks are presented in the concluding section. In the Appendix, the presence of 
the multiplicative anomaly in a very large class of functional determinant regularizations is 
proven. 



II. CHEMICAL POTENTIAL IN THE NON-INTERACTING CASE REVISITED 

We shall treat here the non-interacting bosonic model at finite temperature with a non 
vanishing chemical potential. This situation has been extensively investigated in the past. 
We will reconsider it here making rigorous use of the zeta-function regularization procedure, 
what will allow us to clarify some subtle points overlooked in previous studies. 

Again, it is convenient to work in 5*^ x R^, with D = N + 1 in order to exploit the role 
of the dimension of the space-time. As discussed in the Introduction, here one has to deal 
with the functional determinant of the product of the two operators 

L± = Lr+[VLN±efif , (12) 
while a second factorization is obtained in terms of the operators 




with I/7V = — A/vT + m^. 

Let us start with the latter one, since it renders calculations easier. The zeta function 
can be defined, for a sufficiently large real part of s, by means of the Mellin transform of 
the related heat operator trace, namely 
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1 /■"'-' 

as\Ki) = T,r^ dt f-'Ti exp {-tK±) , (14) 
i [s) Jo 

where Tr exp —tK± is given by 

Tr e"*^± = Tr e"*^^ ^ g-tKiieM)^ _ ^^5) 

n 

For notational simplicity, we write 

Tr e-*^- = ^ e-*^- = -^e"*™' , (16) 
j (47rt) 2 

with Aj the eigenvalues of the operator L^, namely k'^ + m^. Thus, via Mellin transform, 
one has the standard result 

/-/■ \T \ r(2; — y) , , 

(47r)- r(2:) 
The Poisson-Jacobi resummation formula gives 

^ g-t{.„±iM)^ = f 1 + 2 f; cosh(ne/i/?)e-'^^ . (18) 

(Ant) 2 \ J 

As a result, even though ^ K^, one gets 

Tre-*-^+ = Tre~*-^- 

3 ( ^ 2 2 \ 

= Tre-*^-^^-|^ h + 2 cosh(ne/i/3)e-'^ j , (19) 

and C,{s\Kj^ = ({s\K_). Therefore, making use of Eq. (0), Eq. (p!5D, Eq. ([T8| ) and of the 
well known integral representation of the MacDonald function Ks{x), 

K,{Vbx) = -( r e-'^-'^t-'-'dt , (20) 



2 \2Vb, 



one has 



2 ' +^ ' (47r)T r(s) 

+ E cosh(ne/i/3)(n/3)^-2 ^ A, * iT, i (n/^v^,,) . (21) 

The second term in this equation represents the statistical sum contribution, which depends 
non-trivially on the temperature and chemical potential and can be calculated in arbitrary 
dimension. Near s = it can be written as sS{f3,fi) + O(s^), with 

S{p, /i) = - 5] [in (l - e-^(^»+'='^)) + In (l - e"^^^'-^'^))] . (22) 
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Taking Eq. (|T6D into account, this contribution can also be rewritten as follows 



(4vr) 



cosh{ne fj.p){nP) 2 (2m) 2 Kd^uP' 



m] 



(23) 



n=l 



It should be noted that this statistical sum contribution is a series involving MacDonald 
functions. This series is always convergent in any D > 1 space-time dimension and for any 
f3 and e|/i| < m, even in the critical limit e|/ic| = m, m being the lowest eigenvalue in the 
spectrum of Ljv . This follows from the asymptotics of K^{z) for large z: 



(24) 











1 - 


HO(i) 









Coming back to Eq. (21), the first term (the vacuum contribution) has to be considered, 
as usual, for D odd and D even separately. For D odd, we have 



C(o|K+) + c(o|K_) = o 



and 



1 icm^) + c'(oiK-)) = -^n-^)m^ + s{p, /i) . 

2 (An)- 2 

For D even, we shall restrict ourselves to the two cases D = 2 and D = 4. We obtain 



(25) 



(26) 



C(0|ir+)lnM2 + C'(0|i^+) 



An 
' 327r2 



m 



Let us now consider the first factorization. 

Tr e~*"^* = Tr g~*(v^Jv='='=A')^ ^ e 

n 

with 

Tr e~*(^^=^^'')^ = ^ g-*(v^»±eM)2 _ 

i 

Again, Poisson resummation gives 

o / 00 

Tr e~*"^=^ = Tr e~**-^^^^^^^ 
and, as a consequence, 

as\L±) 



(47rt)i/2 



1 + 2 ^ e 4t 

n=l 



(27) 
(28) 

(29) 

(30) 

(31) 



- 2s^ln(l-e 



+ 0(.^ 



(32) 
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Now, for e|//| < m, the binomial theorem yields 



C{z\{VLn + eiif) + C{z\{VLn - eiif) = 2C(^|L 



N 



+ 2 E S^(e/^)^^C(^ + , (33) 



- (2r)!r(2z) 



and this leads to 
1 



{as\L+)+as\L_)) 



(47r) 



r(.) 



(34) 



r(2s + 2r- l)r(g- I) /e^\2r 
^1 (2r)!r(s)r(s- | + r)r(2s- 1) V m J 
- sS{(5,n) + 0{s'') . 



+ E 



r(. + r-^) 



As a result, for D odd, taking the derivative with respect to s and the limit s — > 0, the two 
factorizations give 



lnZ,,,(L+,L_) = ^Yf^V{-^)+S{(5,p.) - \an{L^.L_) , 



(35) 



and 



\nZp,,{K^,K_) = ^^r(-^) + 5(/?,/.) - ^a^(X+,ii'_) , 
(47r)2" 2 ^ 



(36) 



respectively. It should be noted that one gets In Zp^^{L^, L_) = In Zp^^{K^, K_) and thus 
the standard textbook result as soon as one is able to prove that the two multiplicative 
anomalies are vanishing for D odd. We can anticipate that this indeed happens. 

For D even, the situation is different because, within the first factorization, the vacuum 
sector depends explicitly on the chemical potential /i. In fact, for example, one has for D = 2 



I {C'(0|L+) + C'(0|L_) + [C(0|L+) + C(0|L+)] InM^} = ^ 



2/1 -,\ 



+ ^eV + ^(/5,/^), 



(37) 



and for D = 4 

I {C'(0|L+) + C'(0|L_) + [C(0|L+) + C(0|L+)] InM^} 



m 

327r2 



m 



87r2' 



4 4 

^ A'" 2 2 2 

— e n m 



(38) 



As a result, for D = 2 the two factorizations give 



\nZp^^{L+,L_) 



(3Vim 



m 



Vif3 



(39) 
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and 



respectively, while ior D — A one has 



327r2 



m 



m 



+ 



87r2 



e /J, m 



1 



+ S{P,pL)--a4(L+,L^), 



and 



lnZ^,^(i^+,i^_) 



m 



^'(1^^-3/2) 



+ 5(/3,/i)-ia4(ir+,i^_). 



327r2 

In general for D even, one gets 

lnZ^,^(L+,L_) = -/?l^^£^(m,M) +/?V^^£^(m,/x) + - ^az,(L+,L_) , 

and 



lnZ;3,,.(i^+,i^-) = -pVN£v{m,M) + SiP,i^) - -aB(X+,X_) , 
where Sv{m, M) is the naive vacuum energy density, D = 2Q and 



Q 



-2r 



r=l 



Here the cg^r are computable coefficients, given by 

r(2r-l) 2^-^ (-!)«- 



r(2r + l)7rf r(r-|)(g-r)!- 

In the cases D = 2, D = 4 and D = 6 one obtains 



82(771, n) = T^eV^ , S^im, /i) = --^ 



2 2 I 2 e^^^ 
e n [m — 



£6{m,ii) 



1 



167r3 



eV — »^ H ^- — 

^ I 4 6 45 



III. THE INTERACTING CASE IN THE ONE-LOOP APPROXIMATION 



In this section, we will review the interacting charged boson model and we will com- 
pute, within the one-loop approximation, the corresponding operators appearing in the two 
possible factorizations of the partition function. 

The Euclidean action for a self-interacting charged field, again described by two real 
components (pi, is 



S 



dx 



D 



(f), [-Ad + m^) (pi + Ad(0^)^ 



(48) 



where (f^ = (j)k(l)k is 0(2) invariant and Ad = is a dimensionless coupling constant. In the 
important case D = A one has 



S 



dx 



(49) 



Finite temperature effects with non vanishing chemical potential are accounted for in the 
following one-loop partition function p|,3,p!^,p| 



^?i(r)=0(T+/3) 



where 



^0 = /3Vj 



N 



1 



(50) 



(51) 



$ being the background field, assumed to be constant, and A the Euclidean small distur- 
bances operator, given by 



A 



Lij + 2e/ie,jA/Z^+ -^;^%^($2)fcS<l>i$j- , 



{D-2) 



(52) 



with 



(53) 



Again, one is dealing with a matrix-valued, elliptic, differential, non self-adjoint operator 
acting on scalar fields in x R^. As a consequence, the partition function reads [H] 



In Zp{fi) = — - In det 



A 



ik 



M2 



— - In det 



M2 M2 



where 



L± = Lr + Ln + e^/i^ + hD± 



4 



/i^ + 4eV(^7v + ho) 



(54) 



(55) 



with 
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The other factorization is 

K± = Ln + U - + ± (;^^D - 4eV^^r) ' , (57) 

and, again, we have L+L^ = Kj^K^ and we have found the couples of \l/DOs which en- 
ter in the multiphcative anomaly. The necessity of a general formula for computing the 
multiplicative anomalies is clear. This issue will be discussed in the next section. 



IV. THE WODZICKI RESIDUE AND THE MULTIPLICATIVE ANOMALY 

FORMULA 

For the reader's convenience, we will review in this section the necessary information 



concerning the Wodzicki residue |]T6| (see |0 and the references to Wodzicki quoted therein) 
that will be used in the rest of the paper. Let us consider a D-dimensional smooth compact 
manifold without boundary and a (classical) ^'DO, A, acting on sections of vector 
bundles on Mo- To any classical \E'DO, A, it corresponds a complete symbol A{x,k) = 
^-ikx j^^ikx ^ such that, modulo infinitely smoothing operators, one has 

m{x) ~ ^ j^, dy e^^^-y^'A{x, k)f{y) . (58) 

The complete symbol admits an asymptotic expansion for — oo, given by the series 

A(x,A:) ~^A,_,(x,A;) , (59) 

j 

and the coefficients (their number is infinite) fulfill the homogeneity property Aa-j {x, tk) = 
t°-~^ Aa-j{x, k), for t > 0. The number a is called the order of A. For example, in the case 
of differential operators, the complete symbol can be obtained by the substitution 9^ ik^ 
and it has a finite number of coefficients and the series stops at Ao(x, fc), being a polynomial 
in k of order a. 

Now let us introduce the notion of non-commutative residue of a classical \E'DO A of 
order a. If P is an elliptic operator of order p > a, according to Wodzicki one can construct 
the ^DO P4(m), ordP4(M) = p 

Pa{u) = P + uA, (60) 

and its related zeta-function 

Cis\PA{u)) = Ti (P + uA)-' . (61) 

This zeta-function has a meromorphic analytical continuation, which can be determined by 
the standard method starting from the short t asymptotics of Tr exp —tPA{u). Since Pa^u) 



is a ^fDO of order p, the heat-kernel asymptotics reads Wi\ 
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rY^^-tPA(u) ^ Y,a,{u)t'l^ + J2Pkiu)tHnt . (62) 

j=0 k=l 

It should be noted the presence of logarithmic terms in this asymptotic expansion, absent if 
one is dealing with a differential elliptic operator. Taking the derivative with respect to the 
parameter u, one also gets the meromorphic structure of lim^^o '^Ci^lPAiu)). By definition, 
the non-commutative residue of A is 



res (A) = pRes 



\imJ-C{s\PA{u)) 
du 



(63) 



where Res is the usual Cauchy residue. It is possible to show that this definition is indepen- 
dent of the elliptic operator P and that the trace of the operator AP''^ exists and admits 
a meromorphic continuation to the whole complex plane, with a simple pole at s = 0. Its 
Cauchy residue at s = is proportional to non-commutative residue of A: 

res(A) = pRes ,=oTr {AP~') . (64) 

Strictly related to the latter result is the following one, involving the short-t asymptotic 
expansion 

Tr(Ae-*^) ~5]c,t'^-'-^^^^lnt + 0(tlnt). (65) 

Thus, the Wodzicki residue of the \E'DO A can be read off from the above asymptotic 
expansion, selecting the coefficient proportional to Int. Furthermore, it is possible to show 
that Tes{A) is linear with respect to A and that it possesses the important property of being 
the unique trace on the algebra of the classical ^'DOs, namely, one has ies{AB) = Tes{BA). 

Wodzicki has also obtained a local form of the non-commutative residue, which has the 
fundamental consequence of characterizing it through a scalar density. This density can be 
integrated to yield the Wodzicki residue, namely 

res(A)=/ f A^Dix,k)dk. (66) 

Jmd [ztt)^ J\k\=i 

Here the homogeneous component A_u{x, k) of order —D of the complete symbol appears. 
The above result leads to res (A) = when A is an elliptic differential operator. 

Now let us discuss the multiplicative anomaly formula, due to Wodzicki. A more general 
expression has been derived in [l^. Consider two invertible, self-adjoint, elliptic \E'DOs, A 
and B, on M^. If we assume that they commute, then the following equality (the Wodzicki 



multiplicative formula) holds |14| 



res 



(ln(A''5-"))2 



a{A, B) = \ ^. ^. ^ = a(5. A) , (67) 



2ah{a + h) 

where a > and 6 > are the orders of A and 5, respectively. 
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A sketch of the proof is presented in what follows. Recall that if B is an elliptic operator 
of order b > q, according to Wodzicki [|14|, one has the following property for the non- 
commutative residue related to the \1/D0 Q (Eq. ( |6lD is just a consequence of it): in a 
neighborhood of 2 = 0, it holds 



(68) 



Here C is the Euler-Mascheroni constant and tqIB) a coefficient depending also on B. As 
a consequence 



Tr 



Q(B 



-712 



-B 



-72 z 



(72 - 7i) 
267172 



res(Q) +0(z) 



(69) 



in which and 71, 72 are positive real numbers. 

Then it follows that, if is a \1/D0 of zero order and B a \1/D0 of positive order b, and x 
a positive real number, then, in a neighborhood of s = 0, one has r^"^* = 1 — xslnrj + 0(5^) 
and 



sTr 



In r^r] """(5" 



7is 



-72 s 



(72 -7i)res(lnr/) ^^ (72 - 7i)res[(lnr7)^] ^ ^^^^^ 



71726 



72716 



(70) 



As a consequence. 



lim 

s-,0 



sTr(ln 7777-^^(5-^1" - S-T2") = -^(72 - 71) 



res[(ln7/)^ 
72716 



(71) 



Consider now two invertible, commuting, elliptic, self-adjoint operators A and B on M^, 
with a and 6 being the orders of A and B, respectively. Within the zeta-function definition 
of the determinants, consider the quantity 



F{A,B) 



det{AB) 



^a{A,B) 



(det A)(det5) 
Introduce then the family of \l/DOs 

^(a;) = r/^5t, r] = A^B-'', 

and define the function 

det{A{x)B) 



F{A{x),B) 



(det A(a;))(det5) 



(72) 



(73) 



(74) 



We get 



„ a + b 

det5— 



1 _ det(Ag) 
^^^y'^^-(detA)(deti?) 



(detSi)(det5) 

As a consequence, one is led to deal with the following expression for the anomaly 



F{A,B). (75) 
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a{A{x),B) = In F{A{x),B) 



lim 5, 



Tr {A{x)B)~' - Tr A{x)~' - Tr 



(76) 



This quantity has the properties: a{A{0),B) = and a{A{^),B) = a{A,B). 

The next step is to compute the first derivative of a{A{x), B) with respect to x, the result 
being 



dxaiAix)^B) = hm dgS 
Making use of Eq. (^), one obtains 

d,a{A{x),B) 



Tr (\n 



rjrj 



B b 



X 



a{a + b) 



res[(ln?7)^] . 



(77) 



(78) 



And, finally, performing the integration with respect to x, from to 1/6, one gets Wodzicki's 
formula for the multiplicative anomaly, namely 



a{A,B) = a{B,A) 



res 



{IniA'B-"))^ 



2ab{a + b) 



(79) 



It should be noted that a{A, B) depends on a classical ^DO of zero order. Thus, it is 
independent on the renormalization scale M appearing in the path integral. 

We conclude this section by observing that the Wodzicki formula is also valid for \E^DOs 
formally non self-adjoint, provided they are complex functions of self-adjoint elliptic opera- 
tors. An example is given by K±. 



V. THE MULTIPLICATIVE ANOMALY FOR THE INTERACTING 0(2) MODEL 

In this section we come back to the problem of the computation of the multiplicative 
anomaly in the model considered in Sec. 3. Putting ho = 0, which is proportional to the 
coupling constant, we also get the results valid in the free case. Strictly speaking, the results 
of the last section are valid for a compact manifold, but in the case of IR^, or in the finite 
temperature case ^1 X lR^-\ the diver gence is trivial, being contained in the spatial volume 
factor. 

Since the multiplicative anomaly depends on the regularization of the ultraviolet diver- 
gences (in our approach this is equivalent to the asymptotic behavior of the eingenvalues 
at infinity and to the related divergence of the functional determinant) it follows that, due 
to the fact that the Minkowski space-time is ultrastatic {goo = —1), its dependence on the 
temperature is simply proportional to (3, as we will see. 

First, let us consider the self-adjoint factorization. Now the order of L± is 2. Then, the 
Wodzicki formula gives 

an{L+, LJ) = ^ res \{ln{L+Lz')f] . (80) 

There is no ordering problem because we are dealing with commuting operators. We have 
to construct the complete symbol A{x, k) of the classical \I'DO of zero order [ln(L+Ll^)]^. 
The complete symbol reads 
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A{x, k) 



In ( + + e^^ + ho + ( -^h^ + 4eV(^^ + + ho, 



In ( A;^ + + + - ( ;^^d + 4eV^(^^ + + h 



(81) 



where k"^ = k'^ + k^. Let us denote its asymptotics for large |fc| as 



(82) 



According to Wodzicki, we have to select the component and use the multiplicative 
anomaly formula, namely 



1 PVn 



(83) 



For large fc^ and /c^, due to the k"^ dependence, it follows that the homogeneous components 
with odd indices are vanishing. As a consequence, from Eq. (^) one immediately gets the 
following: for D odd, the multiplicative anomaly vanishes. This result is consistent with the 
general theorem contained in | |15|| . 

For D even, the anomaly is present and the asymptotic expansion, from which one can 
easily read off the even homogeneous components, is 



A{x, k) 



e^u^k"^ 

16eV^^ - 32(m2 + eV + ho) + hi + 16eV(^^ + ho) 



+ 



(A;2)2 
128 e^ij\9f 



k^ 



3 (P)2 

Thus, the first non vanishing even components are 

16eV'^' 



(84) 



A_2{X, k) 

A_i{x, k) 



(P)2 
1 



(A;2)2 
128eV'(^')' 



32(m2 + + ho)-^^ + hi + 16eV(^^ + ^d) 



(85) 



It follows that 



a2(L+,L_) = — e /i 

/TT 



(86) 



and 
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a4(L+,L_) 



87r2 



p2 2 
2 2/^ f 



(87) 



If e = 0, we recover the results of Ref. in particular the absence of the multiplicative 
anomaly for D = 2. In the free case h^, = and we have 



a4{L+,L. 



87r2 



2 2/*^ A'' 2\ 

e {— m ) 



(88) 



Let us consider the other factorization. The complete symbol reads now 
A{x, k)= In I A;^ + - + ho + i^^^D - 4eV^^r 

In ( A;^ + - + ho - ( ^^z) " 4eV^^r 



(89) 



For large fc^ and /c^, due to the fc^ dependence, it follows again that the homogeneous 
components with odd indices are vanishing and, for D odd, the multiplicative anomaly is 
absent as in the other case. 

For D even, the first non vanishing even components are 



y4_4(a;, k) 



16eVg 
(A;2)2 



128 eV(^'^' 

3 {k^y 



(A;2)2 

A straightforward calculation gives 



fc2 



PV, 2 2 



(90) 



and 



p'^ll'^ /)2 

eV(m^ + /^i.-^) + ^ 



(91) 



(92) 



Again, for e — > we get the result of ||17 



Some remarks are in order. According to the results of Sec. 2 and this section, it seems 
quite natural to make the conjecture that, in the free case and for D even, one has 

aD{K+,K^) = -aD{L+,L^) . (93) 

We have proved this conjecture for D = 2 and D = A, with an explicit computation, namely 



a2iK+,K. 



m 2 2 

e u 

27r ^ 



p2 2 
2 2/ 2 e ft , 

e fi [m ^ 



(94) 



15 



Thus, in the free case, we also have 

In Z^{L+,L_)= In Zf,{K^,K_) 



(95) 



On the other hand, assuming Eq. (^) and Eq. (p3| ) one obtains a general expression for the 
multiplicative anomaly, namely 



aD{K+,K_) = -az)(L+,L_) = -pViySDim,^) , 
where SDiiTijfi) is given by Eq. (^S]). For instance, for D = 6, we have 



167r3 



4, ,4 



45 



As a consequence, in the free case, one obtains 



where 



In Zp,^{K^, K^) = -(3VNSv{m, M) + VnS{I3, /.) + -/3 W^(m, ^^) 



(96) 



(97) 



(98) 



(99) 



VI. THE FREE CHARGED BOSONIC MODEL AT FINITE TEMPERATURE 

The presence of this new anomaly term in the grand canonical partition requires a re- 
analysis of the statistical mechanical behaviour. This is better accomplished by the use of 
the effective potential formalism at finite temperature and charge. Part of the material of 
this section has appeared before (see for example, p and the references quoted therein) . We 
shall keep working in arbitrary spacetime dimensions, returning to the physically important 
4D case at the end of this section. 

In order to study the effective action in the presence of a non vanishing mean field, it 
is necessary to introduce the one-loop grand partition function in the presence of external 
(constant) sources Jj and the chemical potential /i. This is done with the addition of the 
coupling term — / dx^{Ji(f)i) in the action. As the action is quadratic in the field degrees of 
freedom, a trivial Gaussian integration gives 



Zp{jj, Ji) = exp 



Ji.Ji 



Zpifi) . 



2{m? — e^iJ?) 

where is the spatial volume. In terms of generating functionals per unit volume 



(100) 



N 



In Zi3{fi,Ji] 



1 



N 



In Z/sifi) 



;ioi) 



we have 
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y (A J.) = V^o(A P) - ^^J%^,^ ■ (102) 

which is nothing but the thermodynamics potential of the non interacting boson gas. The 
effective potential in presence of external sources is the Legendre transform of the generating 
functional per unit volume (that is to say, the thermodynamics free energy density), and 
reads 

F(/5, p, $,) = V{I3, p, J,) + /ip + J,<l>„ (103) 

where p and Ji have to be expressed as functions of the charge density, p, and the mean 
field $j, by solving the defining equations 

JV^^<R>^ = 1 /,.«<^.>, (104) 

ap Vn oJi Vn J 

where Vg is a shortened notation for V{(3, p, Ji). As a consequence, with x = |$|, eliminating 
the external sources, one has 

F{P, p, x) = VoiP, p) - ^^MM + 1(^2 + e'ix')x' , (105) 

p = t; h e px , (106) 

ap 



where p is assumed to be a function p = p(/3, p, x), obtained by solving Eq. ( |106| ) in the 
unknown p. 

The possible equilibrium states correspond to minima of the effective potential with 
respect to x, at fixed temperature and charge density. Making use of Eq. ( |105| ) and Eq. ( |106p , 
one arrives at 

1^ = x{m' - eV^) , = K - eV') - 2xeV^. (107) 

ox OX'' ox 

Requiring the first derivative to vanish gives either: 

(i) The unbroken phase, with solution x = 0. From the expression of the second derivative 
this will be a minimum provided ep < m, and the related free energy density will be 

= min^. F{(3, p, x) = Vo{f3, p) - ^^M^li^ = Vo{f3, p) + pp , 

" (108) 

ap 

(ii) The symmetry breaking solution 

ep = ±m , (109) 

here derived as an extremal property of the effective potential. In this case |$| > 0, and 
assuming ep = m, one has 
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Tfl 

= miiia; F{P, p, x) = VoiP, efi = m) + —p 



2 ^dVo{P,p) 
emx = p + 



dp 



(110) 

eii=m 



In this symmetry breaking phase, one has the phenomenon of Bose-Einstein condensation 
(BEC). From the last equation and the fact that Vo{(3,p) is a monotonically decreasing 
function of /3, it follows that = is non vanishing when j3 > (3c ■, the inverse of the 
critical temperature of the BEC transition, as given by the condition 



5K)(/5c,p) 

P = - 



dp 



(111) 

eii=m 



in agreement with the recent criterium for the relativistic Bose-Einstein condensation pro- 
posed in I^O],^. For (3 < (3c, on the other hand, |$| = and the t/(l)-symmetry is restored. 



Now in the free case we know Vg. In the unbroken phase, a; = 0, e\p\ < m, P < (3c, and 
including the anomaly term, Eq. (^5]), with 2Q = D, one has 



D 



4m- ~ / 1 \ 2 _ _ , _ 1 



Q 



Tp = Ev- E cosh(re/i/3)irz, (rm/5) + /xp + - ^ CQ,,(e/x)^^m2«-2- 



P = 7^ E smhirep(3)Ko (rm/5) - e J] CQ,,r (ep)2^-im2«-2^ . (112) 

(27r) 2 Vpr/ 2 

The last term in both equations accounts for the anomaly. 
In the broken phase, e\p\ = m, (3 > (3c, one obtains 

4mT °° f I \ ^~^ 1 ^ o 

Tf3 = Ey- o E cosh(rm/3)/s:D (rm/3) + - E cg.r"^^'^ + — P , 

(27r) 2 vpry z e 

4emT °° / 1 \ "^"""^ 

/5 = ^ E sinh(rm/3)_ft'c (rm/5) — e E "fcq^rf^^^^ + emx^ . (113) 

(27r) 2- V/?ry 2 

In this case, if we formally take the zero temperature limit (3 oo, we get 

Q 

p ^ — e ^ rcQ^r'^^~^ + emx^ . (114) 

r=l 

Now observe that the anomaly contribution leads to some problems, since it does not have 
a definite sign. For example, for any fixed p = 0, in (|114|) could be made negative when 



the anomaly contribution is negative. We thus arrive at a contradiction, as happens, for 
example, for D = 6. 

At this point, however, we should notice that our discussion has involved "unrenormal- 
ized" quantities only, regularized by using the zeta-function procedure. In particular, the 
charge operator has not been properly defined, and the above result indicates that this op- 
erator was not normal ordered relative to the Minkowski vacuum. Now the charge operator 
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appears in the Hamiltonian multiplied by /i. Therefore any normal ordering ambiguity in 
the charge operator gives rise to an ambiguity in the effective action, which must be a linear 
homogeneous function of /i, namely K/x, with a freely disposable constant K. We have 
therefore the freedom to define the generating functional V{j3, /z, J,) up to the linear term 
K^. This will not change the effective action, the Legendre transform being unaffected by 
linear terms. We then introduce the renormalized generating functional 



V'' = V + ^^K = £y + \£o{m, /i) - - . ^. + /zi^, (115) 



where Snifn, ^) is the anomaly, given by Eq. ( ^Sj) and /i) was defined in Eq. (|23|). So 
we have also a finite renormalization of the charge density: 

p« = -^ + e/.x2 = p-K (116) 

Since the effective potential is the same, the extremals are unchanged and we may rewrite 
the free energy in the symmetry breaking phase as 



P e e / 



= ~^^^\e,=m -K- -Y.rcQ.rr?'' + emx^ . (117) 
On the other hand, in the unbroken symmetric phase we have, 

J'p = Sv- ls{f3, p) + pp« + pK + 1 1] CQ,.(e/i)2'■m2«-2^ 

= —-^^ - K - ej:rcQ,{e^r~'m'<^-'^ . (118) 

The constant K can now be fixed by imposing normalization conditions. By doing this 
we see how the spacetime dimension affects renormalization and how special the case D = 4 
is. We do so just by demanding that, at zero temperature and charge density, the symmetry 
be unbroken, which sounds as a very natural normalization condition. Using ( |117[ ), this fixes 
K to be 

Q 

K = _e^rcQ,V^-^ (119) 

As a result, in the broken phase (|e/i| = m) one obtains 

1 m 1 ^ 

Tp = £v- efi = m) + -p^ + - ^ cq,.(1 - 2r)m^'^, 

p e z ^^-^ 

P^ = ~^^^\e,^m + emx\ (120) 



19 



while, in the symmetric phase, one gets 



r=l ^ r=l 



P 



R 



1 

' (3 d^i 



Q 



2r-1^2Q-2r 



;i2i) 



r=l 



r=l 



Note that in this symmetric phase, one can now take the hmit e ^ to reach the usual 
expression of the free energy density for an uncharged boson gas. For general D = 2Q, 
however, the free energy receives contributions from the anomaly term in both phases. On 
the other hand, with the above charge renormalization, the critical temperature is given 
implicitly by (choosing /i > 0) 



P 



. 



;i22) 



eii=m 



which is the usual condition for the critical temperature [^,^ , and it remains unaffected by 
the anomaly. It is important to recognize that the normalization condition just discussed is 
not the only one that is possible. Another physically natural condition would be to demand 
that at zero coupling, the free energy reduces to the free energy of uncharged bosons. This 
then fixes K to be 



K 



r=l 



2Q-1 



(123) 



Now it is the free energy in the broken phase that remains unaffected, while f3c changes. 

Let us consider now the physically important case D = 4. Then, from the expression 
(ff^) for the anomaly, we have in addition the remarkable identity 



E ^C2,, 



1 ^ 

^ r=l 



1 



247r2 



(124) 



This shows that the two normalization conditions just discussed are actually equivalent and 



fix a unique value K 
we obtain 



247r- 



■. As a result, in the symmetry breaking phase and for D = A, 



m 



Ti3 = 8v - -S{(3,m) H p 

p e 



R 



ex 



1 

m 



p dp 



(125) 



e/j=m / 



which means that the broken phase is totally unaffected by the anomaly. On the other hand, 
in the symmetric phase P < Pc, x = and ep < m, the result is 



Sy-^S{P,p)+pp^-"^J^ 



pem" ^ 1 



2 2 2 -^4 4 

e p m — -e p 



l dS{P,p) 
' p dp 



4 



2e>m^ - -eV^ ) + 



em 
12^ 



(126) 
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We observe that the anomaly term (the term in square brackets in ( |126|) ) is present in the 
free energy above Tq. Moreover, it is again vanishing as e — > 0, and the correct expression 
of the free energy density for the uncharged boson gas is recovered. Although the multi- 
plicative anomaly term itself may be not negligible, as a numerical study of it shows, the 
renormalization used in order to preserve the vacuum structure has rendered the anomaly 
term quite harmless near Tc and below. On the other hand, at ultra relativistic tempera- 
tures, T m, it is negligible in comparison with the Planckian T^-term, since the anomaly 
term increases with mass precisely as (however this holds whenever the charge density 
is much less than T^e~^, in order to neglect the Coulomb interaction). Our conclusion is 
that it gives a relevant correction at intermediate temperatures, say of the order of T ^ m, 
which are relativistic for any known massive elementary particle. 



In this paper some deep results of the zeta-function regularization procedure have been 
employed in order to study rigorously the one-loop effective potential for a fixed charged self- 
interacting scalar field at finite temperature. The chemical potential has been introduced 
according to the path integral approach |]T|J^ making use of method I, as discussed recently 
in 1^. Method I leads to the presence of a multiplicative anomaly term, whose existence had 
been completely overlooked in the literature. Using powerful mathematical expressions, this 
anomaly term can be computed exactly. Its explicit form has been obtained here for D = 2 
and D = 4 (while a simple program yields it for any desired value of D.) We have shown 
it to be vanishing for D odd and also its fundamental importance in getting factorization 
invariance (in the operational sense) of the regularized one-loop effective action. On the 
other hand, assuming factorization invariance, we have obtained, in the free case, a general 
expression for the multiplicative anomaly, valid for any even D. 

The existence of this new contribution has led us to revisit and discuss the non-interacting 
case in detail. In particular, we have reexamined the spontaneous symmetry breaking issue 
and the related relativistic Bose-Einstein condensation phenomenon. A renormalization of 
the charge density has been introduced in arbitrary even dimension and has led to the 
usual expression for the critical temperature. In particular, in the broken phase, we have 
shown that only in the physically important case D = 4 the new contribution due to the 
multiplicative anomaly may be absorbed in the charge renormalization process. However, 
in the symmetric unbroken phase, the multiplicative anomaly term gives a non vanishing 
contribution, which has been overlooked in previous investigations. Although it is non 
leading in the ultra high temperature regime, nevertheless it can give relevant correction at 
intermediate temperatures of order T ~ m, m being the mass of the charged boson. 

In the interacting case, even in at the one-loop approximation, it is quite difficult — 
within the first factorization — to deal with the closed expression for the zeta-functions. In 
fact, for D = 4, we have 



VII. CONCLUDING REMARKS 



as\Li) 



/^r(g-i) 
2v^r(s) 



Cis - 2 1^3 + eV + h±J- + 4eV(^3 + h)) 
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2s Tr In 



1 — exp P 



\ 



\ 



Ls + e2^2 + ± j'^ + 4eV(^3 + h)\ + 0{s^) , (127) 



while, within the second factorization, one has difficulties with the sum over the Matsubara 
frequencies Un- For example, again for D = 4, one should deal with 



T(s) Jo 



dtt 



s-l" 



(47rt)^ 



I] exp 



(128) 



As a consequence — as far as the factorization invariance of the partition function is 
concerned — the relevance of the multiplicative anomaly (evaluated in Sec. 4) appears 
manifest, in the interacting case. The issue will require further investigation. 
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VIII. APPENDIX 

In this Appendix we will show that the multiplicative anomaly is present also in a 
large class of regularizations of functional determinants appearing in the one-loop effective 
action. They are often used in the literature and can be called "generalized proper-time 
regularizations", since the prototype among them was introduced by Schwinger [Q. We 
will show that the presence of the multiplicative anomaly stems from the obvious — but 
crucial — fact that the finite part of all these regularization contains, unavoidably, the zeta- 
function regularization contribution C,'{Q\L). 

If we work in the Euclidean formalism, the class of regularization we shall be dealing 
with is defined by |T^,|T^ 



In det 



[jp)^3e{t)) = -J^ dtt-'g,{t)TTe-'^ , (129) 



where gs{t) is a suitable regularizing function, necessary in order to control, for e > 0, 
the (ultraviolet) divergences for small t in the integral and is the mass renormalization 
parameter which renders t adimensional. L is an elliptic operator. We note that the zeta- 
function regularization and the related Dowker-Critchley one Q belong to this class. They 
correspond to the choices 



9lit) = :rAwr^^ 9m = t'- (130) 
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As a consequence, the first regularization, as is well known, gives a finite result ('{e\L), after 
an analytical continuation in e, when the parameter goes to zero, while the second one gives 
(C being again the Euler-Mascheroni constant) 

Indet (A) = -r(^)C(^l^) = -km) - C'(0|^) + CC(0|L) + 0{e) . (131) 

In this case, besides the divergent contribution, the finite part contains the zeta-function 
regularization result. 

Let us consider another class of proper-time regularizations such that the functions geit) 
admit the Mellin transform 

Us) = / dt f-^geit) . (132) 
Jo 

Three popular examples are the ultraviolet cutoff regularization gl{t) = 6{t — e) 

/ T. \ roo r 1 

Indet (^— ] (gm = - dt r'Tve-'^ , g^s) = --e^ , (133) 
the point-splitting regularization g'^{t) = 

Indet (^— j (gtit)) = - dt r^e-fTre"*^ , gt{s) = r(-.)e^ , (134) 
and Pauli-Villars regularization, which in our formalism may be expressed as 

gl{t) = {l-e-^^^y"\ (135) 



Thus 



Indet (A) ^glit)) = -fdtt-' (l - e-'fy-' Tr e-'^ , (136) 



with 



gf(s) = - + ycJe) , (137) 

where Cj{e) are constants which diverge in the limit e — 0. Making use of the Parceval- 
Mellin identity, we may rewrite Eq. (|129|) in terms of a complex integral involving the Mellin 

transform of Tr c^'aT^, namely C{z\L)M'^'^ and (7e(t), i.e. 

Indet (-^) {gM = "7^ / dz r{z)az\L)M'^U-z) . (138) 

VM^y 2m JRez>D/2 

Now the meromorphic properties of C{z\L) are known and we have 
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where Ar are the Seeley-De Witt coefficients (computable) and J{z) is the analytical part 
(normally unknown). We also have C(0|-^) = ^d/2, with Auf2 = when D is odd, since 
we are working in manifolds without boundary. Shifting the vertical contour to the left, 
assuming that the Mellin transform g^{—z) is regular for Kez > and has a simple pole at 
z = given, with residue 1 (see, for example, Eqs. (|133|) and ( |134| )), the residue theorem 
yields 

Indet (^^(t)) = _^'(o|L) +C(0|L) (inM^ + c) 

- J2 ^rUr - y)M^^i - C(0|L) Ine + 0(e) . (140) 

As a consequence, the finite part, coming from the singularity at z = contains ('{0\L) and 
the multiplicative anomaly is also present in this large class of regularizations. The sum 
over r contains the ultraviolet divergences, controlled by e. In the one-loop approximation, 
these divergent terms are removed by the corresponding counterterms. 
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